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Abstract 

For a wide class of multi-user systems, a subset of capacity region which includes the corner points 
and the sum-capacity facet has a special structure known as polymatroid. Multiaccess channels with 
fixed input distributions and multiple-antenna broadcast channels are examples of such systems. Any 
interior point of the sum-capacity facet can be achieved by time-sharing among corner points or by an 
alternative method known as rate-splitting. The main purpose of this paper is to find a point on the 
sum-capacity facet which satisfies a notion of fairness among active users. This problem is addressed 
in two cases: (i) where the complexity of achieving interior points is not feasible, and (ii) where the 
complexity of achieving interior points is feasible. For the first case, the corner point for which the 
minimum rate of the active users is maximized (max-min corner point) is desired for signaling. A 
simple greedy algorithm is introduced to find the optimum max-min corner point. For the second case, 
the polymatroid properties are exploited to locate a rate-vector on the sum-capacity facet which is 
optimally fair in the sense that the minimum rate among all users is maximized (max-min rate). In 
the case that the rate of some users can not increase further (attain the max-min value), the algorithm 
recursively maximizes the minimum rate among the rest of the users. It is shown that the problems 
of deriving the time-sharing coefficients or rate-spitting scheme can be solved by decomposing the 
problem to some lower-dimensional subproblems. In addition, a fast algorithm to compute the time- 
sharing coefficients to attain a general point on the sum-capacity facet is proposed. 

Index Terms 

Polymatroid Structure, Multiuser Systems, Multiaccess Channels, Broadcast Channels, Fairness, 
Successive Decoding, Time-Sharing, Rate-Splitting. 

I. INTRODUCTION 

In the multi-user scenarios, multiple transmitters/receivers share a common communication 
medium, and therefore, there is an inherent competition in accessing the channel. Information 
theoretic results for such systems imply that in order to achieve a high spectral efficiency, the 
users with stronger channel should have a higher portion of the resources. The drawback to this 
is the loss of the fairness among the users. Providing fairness, while achieving high-spectral 
efficiency, is thus a challenging problem. 

A lot of research has addressed this problem and suggested different criteria to design a fair 
system. One of the first criteria is known as max-min measure. In this method, the main effort 
is to maximize the minimum rate of the users, by giving the highest priority to the user with the 
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worst channel. In other words, this method penalizes the users with better channel and sacrifices 
overall efficiency. 

By relaxing the strict condition on fairness, the spectral efficiency can be increased. By 
compromising between fairness and throughput, proportional fairness is proposed in [1]. Based 
on this criterion, the rates of users with a stronger channel can be increased with the cost of 
decreasing the rates of users with a weaker channel. Any change in the rates is acceptable if the 
total proportional increase in the rates of some users is larger than the total proportional decrease 
in the rates of the rest. In fact, by relaxing the strict condition on fairness, the spectral efficiency 
increases. In [2], a criterion based on Nash Bargaining solution in the context of Game Theory 
is proposed. This method generalizes the proportional fairness and increases the efficiency of 
the system. 

All of the aforementioned methods deal with a general multi-user system. However, for a 
wide class of multi-user systems, the capacity region has a special structure that we can exploit 
to provide fairness. Particularly in some multiuser systems, the boundary of the capacity region 
includes a facet on which the sum-rate is maximum (sum-capacity facet). In such systems, one 
can benefit from the available degrees of freedom, and determine the fairest rate-vector on the 
sum-capacity facet. 

As a special case, we consider a class of multi-user systems, in which the whole or a subset of 
the capacity region which includes the comer points and the sum-capacity facet forms a structure 
known as polymatroid. For this class of multi-user systems, the sum-capacity facet has a! corner 
points, where a is the number of users with non-zero power (active users). The sum-capacity 
facet is the convex hull of these comer points. This means that the interior points of the sum- 
capacity facet can be attained by time-sharing among such corner points. As an example of such 
systems, it is shown that the capacity region of multiaccess channels (MAC) with fixed and 
independent input distributions forms a polymatroid [3]. In MAC, the sum-capacity is achieved 
by successive decoding. Applying different orders for the users in successive decoding results in 
different rate-vectors, all with the sum-rate equal to the sum-capacity. The resulting rate-vectors 
correspond to the comer points of the sum-capacity facet. Any point in the convex hull of these 
corner points is on the boundary. In [4], it is proven that the Marton inner bound (see [5]) for 
capacity region of the broadcast channel under fixed joint probability of the auxiliary and input 
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variables, with some conditions, has a polymatroid structure'. As another example, we will show 
that a subset of the capacity region for multiple-input multiple-output (MIMO) broadcast channel 
which includes the corner points forms a polymatroid. 

In [3], the optimal dynamic power allocation strategy for time-varying single- antenna multiple- 
access channel is established. To this end, the polymatroid properties of the capacity region for 
time-invariant multiple-access channel with fixed input distributions have been exploited. In [6], 
the polymatroid properties have been used to find a fair power allocation strategy. This problem 
is formulated by representing a point on the face of the contra-polymatroid (see [3], [7]) as a 
convex combination of its extreme points. 

This article aims at finding a point on the sum-capacity facet which satisfies a notion of 
fairness among active users by exploiting the properties of polymatroids. In order to provide 
fairness, the minimum rate among all users is maximized (max-min rate). In the case that the 
rate of some users can not increase further (attain the max-min value), the algorithm recursively 
maximizes the minimum rate among the rest of the users. Since this rate-vector is in the face 
of the polymatroid, it can be achieved by time sharing among the comer points. It is shown 
that the problem of deriving the time-sharing coefficients to attain this point can be decomposed 
to some lower-dimensional subproblems. An alternative approach to attain an interior point for 
multiple access channels is rate splitting [8], [9]. This method is based on splitting all input 
sources except one into two parts and treating each spilt input as two virtual inputs (or two 
virtual users). By splitting the sources appropriately and successive decoding of virtual users 
in a suitable order, any point on the sum-capacity facet can be attained [8], [9]. Similar to the 
time-sharing procedure, we show that the problem of rate- splitting can be decomposed to some 
lower dimensional subproblems. 

There are cases that the complexity of achieving interior points is not feasible. This motivates 
us to compute the corner point for which the minimum rate of the active users is maximized 
(max-min comer point). A simple greedy algorithm is introduced to find the max-min comer 
point. 

The rest of the paper is organized as follows. In Section |lll the stmcture of the polymatroid 

'Throughout the paper, we deal with the systems where the underlying capacity region or a its subset which included sum- 
capacity facet forms a polymatroid. Apparently, the proposed method can be applied over any achievable region which has the 
similar geometrical structure. In this case, the sum-capacity facet is replaced with maximum-sum-rate facet. 
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is presented. In addition, the relationship between the capacity region of some channels and 
the polymatroid structure is described. Section |llll discusses the case in which the optimal fair 
comer point is computed. In Section |lVl the optimal fair rate-vector on the sum-capacity facet 
is computed by exploiting polymatroid structures. In addition, it is shown that the problem 
of deriving the time-sharing coefficients and rate- splitting can be solved by decomposing the 
problem into some lower-dimensional subproblems. 

Notation: All boldface letters indicate vectors (lower case) or matrices (upper case). det(H) 
and denote the determinant and the transpose conjugate of the matrix H, respectively. M ^ 
represents that the matrix M is positive semi-definite. 1„ represents an n dimensional vector 
with all entries equal to one. is a set of integers E = {I,-- - J-El}, where \E\ denotes the 
cardinality of the set E. The set function / : 2^ — > IZ^ is a mapping from all subsets of E 
(there are a total of 21^' subsets) to the positive real numbers. A permutation of the set E is 
denoted by tt and < i <\E\, represents the element of the set E located in the i*'^ position 

after the permutation. For an a-dimensional vector x = {xi, X2, . . . , Xa] G W and S C E, x{S) 
denotes J2ies ^i- Also, for a set of positive semi-definite matrices Dj, D(S') represents Xlig5 

II. Preliminaries 

A. Polymatroid Structure 

Definition [10, Ch. 18]: Let E = {1, 2, . . . , a} and / : 2^ — > TZ+ be a set function. The 
polyhedron 

B{f,E) = {(xi, ...,xa): x(5) < /(5), V5 C Vx, > 0} (1) 

is a polymatroid, if the set function / satisfies 

{normalized) /(0) = (2) 

{increasing) f{S) < f{T) if S CT (3) 

{submodular) f{S) + f{T)>f{SnT) + f{SUT) (4) 

Any function / that satisfies the above properties is termed as rank function. Note that ([T]) 
imposes 21^' constraints on any given vector {xi, . . . , Xa) G B{f, E). 
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Corresponding to each permutation vr of the set E, the polymatroid B{f, E) has a comer point 
v(7r) G 1Z\ which is equal to: 

^7r(i)(71") = < (5) 

/(Wl),...,7r(z)}) 

-/({7r(l),...,7r(z-l)}) z = 2,...,a 

Consequently, the polymatroid B{f, E) has a! comer points corresponding to different permu- 
tations of the set E. All the comer points are on the facet x(£') = f{E). In addition, any point 
in the polymatroid on the facet x(i?) = f{E) is in the convex hull of these corner points. The 
hyperplane x(£') = f{E) is called as dominant face, or simply face of the polymatroid. In this 
paper, we use the term sum-capacity facet to denote the face of the polymatroid. 

B. Capacity Region and Polymatroid Structure 

For a wide class of multi-user systems, the whole or a subset of the capacity region forms a 
polymatroid structure. As the first example, consider a multiaccess system with a users, where 
the distribution of inputs are independent and equal to p{xi), . . . ,p{xm)- Then, the capacity 
region of such a system is characterized by [11], [12] 

{r e n\\v{S) <I{y; {x,,t G S}\{xi,i G 5^}) VS C , (6) 

where y is the received signal, r represents rate vector, I denotes the mutual information, and 
5*^ is equal to E — S. It has been shown that the above polyhedron forms a polymatroid [3]. 

As the second example, we consider the capacity region of a multiple-antenna broadcast 
system. In the sequel, we show that a subset of the capacity region which includes the corner 
points and sum-capacity facet forms a polymatroid. 

Consider a MIMO Broadcast Channel (MIMO-BC) with M transmit antennas and K users, 
where the A;* user is equipped with receive antennas. In a flat fading environment, the 
baseband model of this system is given by 

yfc = HfcS + Wfc, l<k<K, (7) 

where Hk G C^^y-^'^ denotes the channel matrix from the base station to user k, s E c^'^xi 
represents the transmitted vector, and G C^*^^^ signifies the received vector by user k. The 
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vector Wfc G C^'=^^ is a white Gaussian noise with zero-mean and identity-matrix covariance. 
Consider an order of the users (7r(l), 7r(2), . . . , 7r(Ji )). By assuming that user 7r(i) knows the 
codewords selected for the users 7r(j), j = 1, . . . , z — 1, the interference of the users j = 
1, . . . — over user 7i{i) can be effectively canceled based on dirty-paper-coding theorem [13]. 
Therefore, the rate of user 'ir{i),i = 1, . . . , K, is equal to 

det (ljVfe,7Vfe + H^(i) (j2j>i QttO)) 

r-K{i) = log J J r 7 Y' 

det (In„n, + H^(i) [J2j>i QttU) ) j 
where Qwij) is the covariance of the signal vector to user 7r(j). The capacity region is character- 
ized as the convex hull of the union of such rate- vectors over all permutations (7r(l), 7r(2), . . . , iriK)) 
and over all positive semi-definite covariance matrices Qi, i = 1, . . . , K such that Tr ^Xlilfc Q« j — 
Pt, where Pt denotes the total transmit power [14]. In [15]-[17], a duality between the MIMO- 
BC and the MIMO-MAC is established. In the dual MIMO-MAC, the channel between user 
k and the base station is H|. and the covariance of the power allocated to user k is P^. The 
relationship between and Q^, k = 1, . . . , K, has been derived [16]. The duality is used to 
characterize the sum-capacity of the MIMO-BC as follows 

rsc = max log det 1m,m + ^ H|.PfcHfc . 



Pi, ,Px 
K 



k=l 



s.t. ^Tr(Pfc) <Pt, 



k=l 

Pfc h (9) 

The above optimization problem determines the power allocated to each user in the dual MIMO- 
MAC, and consequently, the power of each user in the MIMO-BC. Note that only a subset of 
users is active and the power allocated to the rest is zero. Equation ^ determines the so-called 
sum-capacity facet. If the cardinality of the set of active users is a, i.e. E = {1, ■ ■ ■ ,a}, the 
sum-capacity facet has a\ corner points corresponding to different permutations of the active 
users. Note that the rates of the non-active users remain zero regardless of the permutation. The 
corner point corresponding to a permutation can be computed using (HJ). Assuming the active 
users are indexed by z = 1, . . . , a, we define 

Di = HlP*H,,2 = l,...,a, (10) 
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where P*, i = 1, . . . ,a, correspond to optimizing matrices in It is shown that the corner 
point in (EJl can be reformulated as [16] 

det (Wf + EjXiD^O)) 

r^(j)=log J r-' « = !,•••, a, (11) 

det (^lAf,A/ + Ej<iD^{i)j 

which is the comer point of the dual MAC. 

Regarding the polymatroid structure of the multiaccess channels and considering the duality 

of the MIMO-MAC and MIMO-BC, we can observe the polymatroid structure of a subset of 

MIMO-BC capacity region which includes the sum-capacity facet. However, to provide a better 

insight about the problem, we introduce a special polymatroid and establish its relationship with 

the capacity region of the MIMO-BC. For a set of positive semi-definite matrices D^, we define 

the set function g as, 

^(^) = logdet(I + D(^)) for S C E. (12) 



Lemma 1 Given g{S) defined in f l72|) . the polyhedron B{g,E) defined as follows is a polyma- 
troid. 

Big, E) = {(xi, . . . , X,) e : x(^) < g{S), V5 C E}. (13) 

Proof: Clearly, g{%) = 0. Assume B ^ and C ^ are two Hermitian matrices. If 
B - C ^ 0, then det(B) > det(C) [14, Proposition 1.2]. Furthermore, if A ^ 0, then [14, 
Proposition 1.3] 

det(A + B + C) ^ det(B + C) 

< -. (14) 

det(A + B) - det(B) ^ ^ 

Using above properties, it is straight-forward to prove © and © for the set function g(.). ■ 

In the set function g{S), define Dj as defined in (flOb . It is easy to verify that the polymatroid 

B{g, E) is a subset of the capacity region of the MIMO-BC. The hyperplane x(_E') = g{E) and 

its comer points (fTTT) are the same as the sum-capacity facet and its comer points. Due to this 

property, we focus on the polymatroid B{g, E) (see Fig. [T]). 
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Fig. 1 

Capacity region of the MIMO-BC and its corner Points. The region OABCD is a polymatroid. 

The line EC is the sum-capacity facet. 



III. The Fairest Corner Point 

As mentioned, in some cases, the complexity of computing and implementing an appropriate 
time-sharing or rate- splitting algorithm is not feasible. This motivates us to compute the comer 
point for which the minimum rate of the active users is maximized (max-min corner point). 
In the following, we present a simple greedy algorithm to find the max-min corner point of a 
general polymatroid B{f,E). 

Algorithm I 

1) Set a = a, S = 0. 

2) Set 7r*(a) as 

7r*(a) = arg min f (E - S - {z}) . (15) 

3) If a > 1, then S < — 5" U {TT*{a)}, a < — a — 1, and go to Step 2; otherwise stop. 
The following theorem proves the optimality of the above algorithm. 

Theorem 1 Let the vector v(7r*) be the corner point of the polymatroid B{f, E) corresponding 
to the permutation n* = (vr*(l), . . . , n*{a)). For any other permutation n = (7r(l), . . . , 7r(a)), 

minu^*(j)(7r*) > mintV(i)(vr). (16) 
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Proof: Assume that in the permutation vr*, the user 6 which is located in position / in the 
permutation tt* ( i.e. 6 = 7r*(/)) has the minimum rate 

W7r*(0(7r*) = min?;^*(i)(7r*). (17) 

i 

Let us define two sets: 

. The set of users located before 7r*(/) in tt*: $ = {7r*(l), . . . , 7r*(/ - 1)}. 

. The set of users located after 7r*(/) in tt*: = {7r*(/ + 1), . . . , 7r*(a)}. 
Using Q, we have 

ve{n*) = f{<^U {9}) -/{<!>). (18) 

In the following, we consider different scenarios which generate new permutations and prove 
that in all cases, (fT6t is valid. 

Case 1. Permutation in $ and ^: By considering (fTSb . it is apparent that any permutation of 
the users in $ and \E' does not change the rate of the user 7r*(/) (see Fig. |2l). 





e 




i J 






e 





Fig. 2 

Case 1. Permutation in $ and 



Case 2. Moving a set of users from \E' to the set Assume a set T of users, T C \E', is 
moved from ^ to the set $ to generate a new permutation tt (see Fig. |3l). The rate of the user 
9 in the new permutation is equal to: 

vei-K) = /($ U T U {9}) - /($ U T). (19) 

From we can show that 

/($ U {9}) + /($ U T) > /($ U T U {9}) + /($). (20) 
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TT* = 


$ 


e 














TT = 




e 


$ -T 



Fig. 3 

Case 2. Moving a set of users from ^ to the set $. 



Using dTSt . (fT9b . and (l20b . we conclude that t'6i(7r) < t;6i(7r*), and therefore, minj f7r(i)(^) ^ 
mini ^'7r*(i)(^*)- 

Case 3. Moving one or more users from the set $ to the set \1/ ( with or without moving some 
users from the set ^ to the set Assume that one or more users move from $ to ^ (with 
or without moving some users from the set $ to the set $) to generate the new permutation vr. 
As depicted in Fig. HI assume that the user u is positioned last in the permutation tt among the 
users moved from $ to (user 7r(l) is positioned first and user 7r(a) is positioned last in the 
permutation tt). 



TT = 














y ^ " 
^ y 
y 

y 

y 




TT = 




e 




V 





\ / 



Fig. 4 

Case 3. Moving one or more users from the set $ to the set * (with or without moving some users from 

the set * to the set $). 



Let n be the set of users located before the user u in the permutation it. Using ©, we have, 

vM = fi^^{^})-fi^)- (21) 
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It is clear that, 

{0} U $ - {u} C n. (22) 
Using © with S = ^ U {6} and T = Q, and regarding (ES), we have, 

f{n u {z.}) - f{n) < /($ u {^}) - /($ u {e} - {u}). (23) 

On the Other hand, the user u is in the set $ in permutation vr*. It means that in Step 2 of 
the algorithm, this user has been compared with other users in the set $ U {9} to be located 
in the position /, but the user 9 has been chosen for the position, i.e. / ($ U {9} — {9}) < 
/($ U {9} - {u}), therefore, 

f{<^)<f{<^U{9}-{u}). (24) 

Using (fTSt . (I2T1) . (l23l) . and (l24l) . we conclude that t^,y(7r) < t^6i(7r*), and therefore, we have 
minj f ^(j) (tt) < min^ tv.(j)(7r*). Note that the permutation of users located before (or after) the 
user v in the permutation tt does not increase Vui'n'). ■ 

Remark: For multiple access channels, the above algorithm suggests that to attain the fairest 
comer point with successive decoding, at each step, one should decode the strongest user (the 
user with the highest rate, while the signals of the remaining users are considered as interference). 
Note that in MAC, the comer point corresponding to the specific permutation tt is obtained by 
the successive decoding in the reverse order of the permutation. 

It is worth mentioning that by using a similar algorithm, one can find the corner point for 
which the maximum rate is minimum. The algorithm is as follows: 

Algorithm II 

1) Set a = 1, S = 0. 

2) Set 71* (a) as 

IT* (a) = arg max / (5 + {z}) . (25) 

3) If a < a, then S < — S U {7r*{a)}, a < — a + 1, and go to Step 2; otherwise stop. 

The optimality of the above algorithm can be proven by a similar method as used to prove 
Theorem [H 
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IV. Optimal Rate-Vector on the Sum-Capacity Facet 

A. Max-Min Operation over a Polymatroid 

In the following, the polymatroid properties are exploited to locate an optimal fair point on 
the sum-capacity facet. For an optimal fair point, the minimum rate among all the users should 

be maximized (max-min rate). For a sum-capacity of rsc, a fair rate allocation would ideally 

rsc 

achieve an equal rate of — for the a active users. Although this rate- vector is feasible for some 

a 

special cases (see Fig. |5l), it is not attainable in the general case (see Fig. |6l). The maximum 
possible value for the minimum entry of a vector x, where x E B{f, E), can be computed using 
the following lemma. 




Fig. 5 Fig. 6 

All-Equal Rate-Vector Is on the Sum-Capacity All-Equal Rate-Vector Is NOT on the 

Facet Sum-Capacity Facet 



Lemma 2 In the polymatroid B{f, E), define 

5 = max min Xj. 

s.t. Gi3(/,E). (26) 

Then, 

= mm -777]-. (27) 

ScE,Sj^0 \S\ 
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Proof: Consider x G E), and let a = miniXi. Therefore, 

yS C E, a\S\ < x(5). (28) 

Noting VS* C E, x(S') < f{S) and using the above inequality, we have 

V5 C E, a\S\ < f{S). (29) 

Consequently, a < minscE, 5^0 Therefore, minscE, s^(l> ^j§^ provides an upper bound on 
minjXj. By selecting x = 51a E B{f,E), where 5 = min^^^^;^ 5.^0^^, the upper bound is 
achieved, and the proof is completed. ■ 

In minimization dZTt . if the minimizer is not the set E, then 6 (the optimal max-min value) is 
less than ^ ( rsc = f{E) is the sum-capacity), and therefore, the ideal fairness is not feasible. 
For example, in the polymatroid depicted in Fig |6l the minimizing set in (l27b is the set {3}, 
and therefore 6 = /({3}). 

In the following, a recursive algorithm is proposed to locate a rate vector x* on the sum- 
capacity facet which not only attains the optimal max-min value 6, but also provides fairness 
among the users which have the rates higher than 6. The proposed algorithm partitions the set 
of active users into t + 1 disjoint subsets, S^^\ S^^\ such that in the z'th subset the rate of 
all users is equal to m^''\i = 0, ■ ■ ■ ,t, where S = m^^^ < m^^^ < ■ ■ ■ < m^*\ Starting from 
m^'^\ the algorithm maximizes m^^\ i = 1, - ■ ■ ,t, given that m^-'^'s, j = 0, ■ ■ ■ , i — 1, are already 
at their maximum possible values. To simplify this procedure, we establish a chain of nested 
polymatroids, E^")), a = 0, . . . , t, where 

eW) c B{f'-^\E^'-^^) c . . . C E(°)) = B{f, E). (30) 

In this algorithm, we use the result of the following lemma. 

Lemma 3 Let _E = {1, . . . , a} and A G E, A ^ E. If the set function f : 2^ — > 71+ is a rank 
function, then h : 2^~^ — > 71+, defined as 



h{S) = f{SuA)-f{A), ScE-A, (31) 



is a rank function. 



Proof: By direct verification. 
Using the following algorithm, one can compute the rate- vector x*. 
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Algorithm III 

1) Initialize the iteration index a = 0, E^^^ = E, and f^^^ = f. 

2) Find m'^"\ where 

(a) ■ f^^'H^) 

' = mm — --- — . (32) 

ScE('^\Sj^(l> \o\ 

Set S^"^ equal to the optimizing subset. 

3) For all i E 5^"), set x* = m^"). 

4) Define the polymatroid where 

and C 

= u 5^°)) - (34) 

5) If £'("+1) ^ 0, set « < — a + 1 and move to step 2, otherwise stop. 

This algorithm computes the optimization sets S^°'\ a = 0, - ■ ■ ,t and their corresponding 
m("), where E = Uj^g S^^'^ and x* G {m(°), ■ ■ ■ , m^}, z = 1, ■ ■ ■ , a. 

To provide better insight about the algorithm, let us apply it over the polymatroids depicted 
in figures |5l and |6l For the polymatroid in Fig. |5l the algorithm results in x* = (m^°\ m^^^ , m^^^ ) 
where m*^°) = liiMiM por the polymatroid shown in Fig |6l the resulting point is x* = 
(m(i),m(i),mW), where = ^ and m^) = /^MlM = /({i,2.3»-/({3}) ^^^^ 

In the following, we prove some properties of the vector x*. 

Theorem 2 Assume that the algorithm III is applied over the polymatroid B{f, E), then 
(I) X* G B{f,E) and is located on the sum-capacity facet x(_E) = f{E). 
(II) The minimum entry of the vector x* attains the optimum value determined by Lemma |2] 
and 

5 = m(°) < m(^) < ■ ■■ < mW. (35) 

Proof- 
Fart (I): We show that x* G B{f, E). According to the algorithm, we have m*^°) = rmn.scE,si^% ^f^^ 
where S^^^ is the minimizing set. In addition, x* = m*^°^ for all i G S^'^\ It is straight- 
forward to check that the assigned values for x*,i G S^^\ do not violate the constraints of 
the polymatroid B{f, E), expressed in ([T]). By substituting the assigned values for x^, i G S'^^\ in 
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Fig. 7 

The Fairest Rate Vector x* on the Sum-Rate Facet of the Polymatroid 



the constraints of the polymatroid E), the constraints over the coordinate i, i E E — S^^\ 
are updated as follows: from the definition of the polymatroid, we have a set of constraints on 
x(5'), S <Z E - S^^\ which has the following format: 

VA C 5(°\x(5U A) < /(°)(^U A). (36) 

Since S n A = 0, then x{S U A) = x{S) + x(A). Consequently, from (|36l), we have, 

VA c ,x{S) < f°\S U A) - x( A) . (37) 

Consequently, \/ S C E - S^^\ 

x(5) < min {j<^\S U A) - x(A)}. (38) 
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We claim that min^^smif^^HS U A) - x(A)} is equal to f^^\S U - The proof 

is as follows: 

VA C , f°\S UA)- x(A) (39) 

> /(°)(^U A) (40) 
>;(o)(5U5(o))_^(o)^^(o))_ (41) 

The first inequality relies on the fact that \/A, x(y4) < f^^\A). The second inequality is proven 
by using Q and the fact that A C S^'^^ and S n S^^^ = 0. It is easy to check that the above 
inequalities change to equalities for A = S^^\ 

Regarding the above statements, for the non-allocated entries of x, we have the following set 
of constraints, 

V5 C ^ - x(5) < /(°) ('^ U ) - /(°) ) . (42) 

Let us define E^^^ = E^^) - ^(o), f^^\S) = U S^) - /^(^(o)), \/S C E^^\ By using 

Lemma El the set of constraints (B^ on defines the polymatroid B{f^^\ E^^^), which is 
a subset of B{f,E). Now, we use the same procedure that is applied for B{f^^\E^^^) over 
B{f^^\E^^^), and continue recursively. Therefore, in iteration indexed by a, a = 0, . . . , t, the 
rates of a subset of coordinates are determined such that the constraints of the polymatroid 
S(/(°),E(")) are not violated. Since E^")) c B{f,E), then x* e B{f,E). Direct 

verification proves that x*(E) = f{E). 

Part (II): We must show that the smallest entries of x* is equal to miuscE ^^j§^- According to 
the algorithm, for all i ^ E, we have x* G {m^^\ ...,m'^^^}. Furthermore, m^^^ = minscE "^j^- 

From the algorithm, we have 

(,•) _ f^'Ks^^^) _ . f^^\s) /o-)(go-+i)u^o-)) _ u g(i)) 

Therefore, 

u) /U)(50-+i)u5b))-mO)|Sb)| 
m ^ |5{i+i)| 

where (l46b relies on LHS of (l43t . Consequently, m'^^) < m^^^ < . . . < m*^*) and the proof is 
complete. ■ 
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The remaining issue in Algorithm III is how to compute mmscE,s^$ ^g^- These types of 
problems are known as geometric minimizations. In order to find the minimizer, the smallest 
value of P is desirable such that there is a set S with f(S) = P\S\. For the special case of single 
antenna Gaussian multiaccess channels, computing such /? is very simple. For the general case, 
P can be computed by Dinkelbach's discrete Newton method as follows [18]. 

The algorithm is initialized by setting /? equal to f{E)/\E\, which is an upper bound for 
optimum /5. Then, a minimizer Y of f{S) — P\S\ is calculated, as will be explained later. Since 
f{E) - (3\E\ = 0, then f{Y) - (3\Y\ < 0. If f{Y) - (5\Y\ = 0, the current (3 is optimum. If 
f{Y) — P\Y\ < 0, then we update P = f{Y)/\Y\, which provides an improved upper bound. By 
repeating this operation, the optimal value of P will eventually be calculated [18]. It is shown 
that the number of P visited by the algorithm is at most \E\ [18]. 

Using this approach, the minimization problem ming^e.^^e '^j^ is changed to mmscE,s^0 fi^)" 
P\S\. By direct verification of Q, it is easy to see that f{S) — P\S\ is a submodular function. 
There have been a lot of research on submodular minimization problems [18]-[20]. In [19], 
[20], the first combinatorial polynomial-time algorithms for solving submodular minimization 
problems are developed. These algorithms design a strongly polynomial combinatorial algorithm 
for testing membership in polymatroid polyhedra. 

B. Decomposition of the Time-Sharing Problem 

In the following, we take advantage of the special properties of x* and polymatroids to 
break down the time- sharing problem to some lower dimensional subproblems. In the previous 
subsection, a chain of nested polymatroids B{f^°'\ E^"^), a = 0,...,t, is introduced, where 
c for a = Since C E'^^^ for j = 0,...,t and 

regarding the definition of polymatroid, B{f^^\ S^^^), j = 1, . . . , t, is a polymatroid, which is 
defined on the dimensions S''-'\ According to the proof of Theorem |2l the vector m^^^l^gu)^ E 
B{f^^\S^^'^) is on the hyperplane x(5(^)) = f{S^^^). Let {tt^^^Tj = 1, . . . ,\S^^'^\\} be the set 
of all permutations of the set S^^\ and be the corner point corresponding to the 

permutation 7r[f in the polymatroid B{f'^^\ S^^^). Then, there exist the coefficients < X^^' < 1, 
7j = such that 



|S(J)|! 



m 



(47) 



7,=1 
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where 

E = 1- (48) 

7,=1 

Note that E = IJ*=o Consider a permutation vr-^p as one of the total \S^^^\\ permutations of 
5"*^-^^ for j = 0, ■ • • ,t, then the permutation tt formed by concatenating these permutations, i.e. 
TT = ^vr^t^ ■ ■ ■ , vr^o^ j , is a permutation on the set E. 

Theorem 3 Consider the permutation n = , ■ ■ ■ , n^J ^ of the set E. 
(I) The corner point corresponding to the permutation n in the polymatroid B{f, E) is 

t,.(7r) = np^(7r(^)), for i e S^^\ (49) 

where u(^')(7r^^'^) is the corner point of the polymatroid B{f^^\ S^^^) corresponding to the 
permutation Hyf, and Ui^\TT^l^) denotes the value o/u^-'^(7r^j'') over the dimension i, i E S^^\ 
(II) The vector x* is in the convex hull of the set of corner points corresponding to the following 
set of permutations 

{ (ttW, ■ ■ ■ , vrW) , 1 < 7. < |^(*)|!, . . . , 1 < 70 < \S^'^\l} , (50) 

where the coefficient of the corner point corresponding to the permutation n = (^n^t , ■ ■ ■ , n^J 
is equal to \^^} . . . X^^, i.e. 

|SW|! 

-•=E---E<'---^WW?.---,<:')). (51) 

7t=l 70=1 

Proof: Part (I) From recursive equation we can show that 

For S e ^ - U ^^'^ /^'^(^) = / (^^ {U ^^'^}) " ^ ({[i^^'^}) ■ ^^^^ 

Consider the permutation n = (^n^t , ■ ■ ■ , n^J^ . Set ^ = ^^^^ \S^^''\. By using ^ and (l52t . for 

IS^-^'+^^l, w,r(K)(vr) is equal to 

(vr) = / ({7r(l), . . . , 7r(K)}) - / ({7r(l), . . . , 7r(/s: - 1)}) (53) 



rr 

f^^^ ({7r(e + 1) . . . , 7r(K)}) - f^^> ({7r(e + 1) . . . , vr(/^ - 1)}) . (55) 



.1=0 J / \ K.i=0 
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According to definition of polymatroid and its comer points, the RHS of S55h is the value of 
the corresponding corner point of the polymatroid B{f^^\ S^^^). 
Part (II) Since J2y,!=i ^7? = 1 and by using (l47t and part (I) of the theorem, it is easy to 
verify that the i E S^^\ entry of 

|5(0)|! 

^A^vK;),...,.*?) (56) 

70=1 

is equal to m(°\ Similarly, the entry i, i E S^^\ of 

15(1)1! 15(0)1! 

E^i;'E^rvW? '^5?), (57) 

71=1 70=1 

is equal to m^^\ while the entry i, i E S^^\ remains m^^\ By continuing this procedure, part 
(II) of the algorithm is proven. ■ 

Regarding the above statements, the problem of finding time-sharing coefficients is decom- 
posed to some lower dimensional subproblems. In each sub-problem, the objective is to find 
the coefficients of the time-sharing among the corner points of the polymatroid B{f'^^\ S^^'^), 
j = 0, . . . , t, to attain m^^^l^gu)] - In this part, we present an algorithm which finds the coefficients 
of the time-sharing over the corner points of a general polymatroid B{f, E) to attain a vector x 
located on the face of the polymatroid. 

Algorithm IV 

1) Initialize « = 1, ui = v(7r*) (the fairest corner point obtained by algorithm I). 

2) Solve the linear optimization problem 

maxr 

S-^- Yh=1 fJ'i^i - X > r 

0<fii<l (58) 
Let yu°, z = 1, . . . , a be the optimizing coefficients. 

3) If X = x;r=i i^i^i^ Stop. 

4) a < — a + 1. Set e = x — ^"^^ /i"uj and determine the permutation tt for which 6,^(1) > 
e7r(2) > • • • > ^Trd^i)- Set Uq, = v(7r) and move to step 2. 

The idea behind the algorithm is as follows. In each step, the time-sharing among some comer 
points is performed. If the resulting vector is equal to x, the answer is obtained; otherwise a 
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permutation vr is determined such that ejr(i) > e^(2) > . . . > g-kHei), where the error vector 
e represents the difference between the vector x and resuhing vector from time-sharing. We 
can compensate the error vector e by including an appropriate comer point in the set of comer 
points participating in time-sharing. Clearly, the best one to be included is the one which has 
the highest possible rate for user 7r(l) and lowest possible rate for user Trdi^l). Apparently, this 
corner point is v(7r), computed by algorithm IV. 

Note that Algorithm IV can be applied over the sub-polymatroids B{f^^\ S^^^), j = 0, . . . ,t, 
to attain m^^^l^g^j)^ or directly applied over the original polymatroid to attain x* . If a and \S^ 
are relatively small numbers, the decomposition method has less complexity, otherwise applying 
Algorithm IV over the original problem is less complex. 

C. Decomposition of Rate-Splitting Approach 

As mentioned, an altemative approach to achieve any rate-vector on the sum-capacity facet 
of MAC is rate splitting [8], [9]. This method is based on splitting all input sources except one 
into two parts, and treating each spilt input as two virtual inputs (or two virtual users). Thus, 
there are at most 2a — 1 virtual users. It is proven that by splitting the sources appropriately and 
successively decoding virtual users in a suitable order, any point on the sum-capacity facet can 
be attained. 

Similar to the time-sharing part, we prove that to attain the rate vector x*, the rate- splitting 
procedure can be decomposed into some lower dimensional subproblems. Consider a MAC, 
where the capacity region is represented by polymatroid E) and the vector x*, derived 
in Algorithm III, is on its face. Assume that the users in the set S^^^ are decoded before the 
set of users in {S^^~^\ S'^^"^\ ^(°)} and after the users in the set {^W, . . . , S^^+^\ 3'^^+^^ 
.Therefore, by similar discussion used in (l36b to (l42b . we conclude that the rate of the users in 
the set S^^^ is characterized by the polymatroid B{f^^\ S^^^), where the rate-vector m^^^l^ga)^ is 
on its face. Regarding the results presented in [8], [9], we can attain the rate-vector m(-5'^l|5.(j) 
by properly splitting the sources of all inputs, except for one, in the set S^^^ to form 2|S''^-'^| — 1 
virtual users and by choosing the proper order of the decoding of the virtual users. Consequently, 
using algorithm V (below), we achieve the rate-vector x* in the original polymatroid. 

Algorithm V 
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1) Apply rate-splitting approach to attain the rate- vector m^^'^l^gij)^ on the face of the poly- 
matroid B{f^\ S^^^), for j = 0, . . . , t. Therefore, for each j, < j < t, at most 2\S^^^- 1 
virtual users are specified with a specific order of decoding. 

2) Starting from j = t, decode the virtual users in the set S^^'^ in the order found in Step 1. 
Set j ^ j — 1. Follow the procedure until j < 0. 

V. Conclusion 

We considered the problem of fairness for a class of systems for which a subset of the capacity 
region forms a polymatroid structure. The main purpose is to find a point on the sum-capacity 
facet which satisfies a notion of fairness among active users. This problem is addressed in cases 
where the complexity of achieving interior points is not feasible, and where the complexity of 
achieving interior points is feasible. For the first case, the corner point for which the minimum 
rate of the active users is maximized (max-min corner point) is desired for signaling. A simple 
greedy algorithm is introduced to find the optimum max-min comer point. For the second case, 
the polymatroid properties are exploited to locate a rate- vector on the sum-capacity facet which is 
optimally fair in the sense that the minimum rate among all users is maximized (max-min rate). 
In the case that the rate of some users can not increase further (attain the max-min value), the 
algorithm recursively maximizes the minimum rate among the rest of the users. It is shown that 
the problems of deriving the time-sharing coefficients and rate- spitting scheme can be solved by 
decomposing the problem to some lower-dimensional subproblems. In addition, a fast algorithm 
to compute the time-sharing coefficients to attain a general point on the sum-capacity facet is 
proposed. 
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